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I. INTRODUCTION 

The pion-nucleon interaction has provided increasingly 
sensitive tests of nuclear theory. One of the challenges yet 
to be met is to understand the polarization observables 
for pion production in pp collisions. This is especially 
interesting near threshold where few partial waves con- 
tribute and where calculations should be more manage- 
able and more conclusive. 

After the initial theoretical work in the fifties by Gell- 
Mann and Watson |jj and Rosenfeld ||, more than a 
decade elapsed before explicit pp — ► ppir and pp — > pnir + 
cross sections for Ss (Inn = 0, l w = 0) transitions 
were predicted by Koltun and Reitan |J in 1966 and 
by Schillaci, Silbar and Young Q in 1969. When the 
small cross sections very close to threshold could finally 
be measured twenty years later H|), ^ turned out that 
these calculations had missed the true cross sections by 
factors up to five. This realization has spurred much new 
theoretical research. 

To date the Julich meson exchange model 0-|ll| has 
yielded the most successful calculations. This model rep- 
resents a much advanced development of the approach 
of Ref. [| and builds on the insights of the 1990's (e.g., 
by Lee and Riska |T|] and many others). It has per- 
mitted detailed calculations beyond l„ = transitions 
and provides analyzing powers and spin correlation coef- 
ficients for the near-threshold region. The Jiilich model 
incorporates all the basic diagrams, realistic final state 
interactions, off-shell effects, contributions from the delta 
resonance and the exchange of heavier mesons. With the 
exception of the heavy meson exchange term there are no 
adjustable parameters. At this time it is the only model 
with predictions that can be compared to our measure- 
ments. However, the Julich model does not account for 
quark degrees of freedom, the potential study of which 
had motivated our experiment initially. 

Ideally, one would interpret the basic pion production 
reactions in a framework compatible with QCD, e.g., cal- 
culations using chiral perturbation theory (\PT). How- 
ever, with one exception the xPP calculations pub- 
lished to date are still restricted to — 0. Moreover, 
for all three pp — > Xn reactions, the \PT cross sections 
remain a factor of two or more below experiment [fl4"|| . 
This shortcoming may be attributable to the difficulties 
of xPT f° r momentum transfers larger than m w . The 
xPT calculations published to date are best viewed as 
work in progress [p.5| . 

Calculations and experiments very close to threshold 
require great care. For pp — > pnir + Ss transitions (l n = 0, 
Ipn = 0) only one amplitude is calculated and the angu- 
lar dependence is trivial. However, the near-threshold 
cross section and its energy dependence are significantly 
modified by "secondary" effects, such as final state inter- 
actions that are particularly important for l pn = 0. 

Measurements very close to threshold can present dif- 
ficulties because the cross sections are small, of the order 



of 1 fib, and change rapidly with energy. The energy of 
the interacting nucleons for reactions very close to thresh- 
old must be precisely known and maintained. At IUCF 
this has been accomplished by the use of a very thin in- 
ternal target and the precise beam energy control of the 
Cooler (storage) Ring. The IUCF Cooler also generates 
low background. 

The earliest studies of pp — > pnn + very close to 
threshold |p|,^6[-^8[ had available a stored IUCF beam 
of < 50/iA They used an unpolarizcd gas jet target and 
measured cross sections and analyzing powers from 293 
MeV (i.e., 0.7 MeV above the tt + production threshold) 
to 330 MeV. These experiments deduced cross sections 
for Ss pion production very close to threshold. As long 
as Ss production of pions strongly dominated, analyzing 
powers also provided information for Sp (l n — 1) admix- 
tures |fl8|] . At 325 MeV and above, higher partial waves 
enter significantly, but the larger cross sections make it 
practical to explore analyzing powers and spin correlation 
coefficients, which allow a much more detailed compar- 
ison of theory and experiment. At the upgraded IUCF 
Cooler Ring, an intense polarized proton beam with a 
large longitudinal component and an efficient window- 
less polarized hydrogen target now permit measurements 
of all spin correlations coefficients for pp — > pmt + . Some 
initial results for transverse spin correlations have been 
reported in [ pj| . 

The goal of the present study is to quantify the grow- 
ing importance of higher partial waves (Sp, Ps, Pp, and, 
potentially, Sd transitions) by measuring analyzing pow- 
ers and spin correlation coefficients as a function of en- 
ergy. These polarization observables are sensitive indi- 
cators of the reaction mechanism and the contributing 
partial waves Jl8],^0 21 and are a powerful tool in deter- 
mining transition amplitudes empirically. In this exper- 
iment we measured A xx + A yy , A xx - A yy , A zz , A xz , 
as well as the polarization observables A y and A z for the 
energy region 325 to 400 MeV. 



II. EXPERIMENT 
A. Experimental considerations 

The Cooler Ring of the Indiana University Cyclotron 
Facility (IUCF) produces protons of energies up to 500 
MeV with polarization of P 0.65, low emittance and 
low background. This permits in-beam experiments of re- 
actions with microbarn cross sections. The improvement 
of beam intensity at IUCF over time now allows the use 
of very thin polarized targets. During the pp — > pmr + ex- 
periment typical intensities of the stored polarized beam 
ranged from 100 — 300/iA 

The apparatus for polarized internal target experi- 
ments (PINTEX) makes use of a windowless target cell 
continuously filled by a polarized atomic Hydrogen beam. 
The measurements cycle through a full set of relative 
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beam and target spin alignments. The technical aspects 
of beam preparation, electronics, and target and detector 
properties have been reported previously in p2[ . As is 
customary, the beam is defined to travel in the positive 
z direction, y is vertical and x completes a right handed 
coordinate system. Below is a brief review of parameters 
pertinent to the data analysis. 




FIG. 1. The PINTEX detector for the experiment: F is a 
thin timing detector. S labels one of the four detectors for 
the elastic pp scattering monitor. WC1 and WC2 are wire 
chambers. E and K are segmented plastic scintillator stacks 
that determine the energy of the charged reaction products. 
V is the charged particle veto detector, and H is the neutron 
hodoscope. 



nal (z) and vertical (y) spin components and a very small 
component in the (x) direction, as shown in Table I. 

Elastic pp scattering was used to measure and monitor 
the three beam polarization components as well as the 
luminosity. Elastic protons were detected with four plas- 
tic scintillators mounted at 9 — 45°, with cf) = ±45° and 
±135°. Coincident protons striking these monitor detec- 
tors (labeled S in Fig. 1) pass through wire chamber 1, so 
the needed tracking information is available. The prod- 
uct PQ of beam polarization (P) and target polarization 
(Q) was deduced from the large known spin correlation 
A xx — A yy in elastic scattering |2f| . A 3D sketch of the 
detector system is shown in Fig.l. 

The reaction pions in this study had lab energies from 
0.1 to 120.5 MeV and were emitted at polar lab angles 
from 0° to 180°. By contrast, the reaction nucleons re- 
main constrained by kinematics to forward angles below 
31.2° and to lab energies from 20.8 to 227.9 MeV. This 
range of angles and energies affects the choice of detec- 
tors that can be employed. If both outgoing nucleons are 
protons as in pp — > ppir one can ignore the pion and use 
a moderate size forward detector to intercept almost all 
ejectiles of interest |^2). This procedure was used for the 
simultaneously measured pp — > ppir reaction |26| ] . The 
corresponding procedure for pp — > pnn + is to mount a 
large area neutron hodoscope behind the proton detec- 
tors and determine the energies of the detected neutrons 
by time of flight. The construction and operation of the 
neutron hodoscope has been described previously in p7| 



In the experiment we used the Madison atomic beam 
target with a storage cell of very low mass |23|,Q. The 
storage cell had a length of 25 cm and a diameter of 1.2 
cm. This open-ended cylindrical cell produces a trian- 
gular shape of the target density distribution with its 
maximum at the center (z=0). It was made of a thin (25 
fim) aluminum foil to keep to a minimum background 
events caused by the beam halo. Teflon coating was used 
to inhibit depolarization of the target atoms. Sets of 
orthogonal holding coils surround the storage cell. The 
coils are used to align the polarized hydrogen atoms in 
the ±x,±y, and ±z directions. Typical target polariza- 
tions were Q = 0.75 and the approximate target density 
was 1.4 x 10 13 atoms/cm 2 . 

The target spin alignment can be changed in less than 
10 ms. During runs the target polarization direction was 
changed every 2 sec and followed the sequence ±x, ±y, 
and ±z. Each data taking cycle had constant beam po- 
larization and was set to last 5-8 minutes, after which 
the remaining beam was discarded. The beam polariza- 
tion was reversed with each new cycle to minimize the 
effect of apparatus asymmetries. In the first phase of the 
experiment (run a) the beam spin directions were alter- 
nated between +y and -y. In the more recent runs (b) 
solenoid spin rotators were used to give the beam spin 
a large longitudinal component. This spin rotation was 
energy dependent and produced roughly equal longitudi- 
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FIG. 2. Raw ejectile time of flight (ADC channels) vs. 
energy deposited in the E detector by pions, protons and 
deuterons. Triggers from two charged particles tracks and 
(no K) were a prerequisite. This spectrum was used for iden- 
tification of protons and pions. 

All detectors are segmented because the energies of the 
coincident reaction particles need to be measured inde- 
pendently. Monte Carlo calculations suggest that eight 
Aij!> segments are sufficient because of the tendency of 
the ejectiles to have significantly different azimuthal an- 
gles. The K detector was needed to obtain the necessary 
stopping power for the more energetic pions and protons. 
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Identification of the charged particles was usually accom- 
plished by their time of flight vs. energy correlation, 
where the start signal was supplied by the F detector 
and the stop signal was provided by the E detector. The 
pion and proton distributions were generally well sepa- 
rated. Fig. 2 shows a typical particle ID spectrum for 
accepted pn + coincidence events. 

The more energetic ejectiles stop in the K detector, and 
superior particle identification is obtained by comparing 
the energies deposited in the K vs. the E detectors, as 
seen in Fig. [| 
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FIG. 3. Particle identification cuts at 375 MeV for ener- 
getic ejectiles based on energy deposited in the K detector as 
function of energy loss in the E detector. Acceptable events 
had to be inside the regions outlined. 

We measured the polarization observables Aij in two 
different ways: 

1) by measuring the pions directly, in coincidence with 
protons (p+7r method), and 

2) by reconstructing pion momenta from the measured 
proton and neutron momenta (p+n method). 

The first method had the advantage of simplicity and 
high count rate, but we cannot measure pions at large 
angles due to the limited detector size. Therefore the spin 
dependent cross section ratios could be compared with 
theoretical spin correlation coefficients only at forward 
angles. The second method is free from this limitation, 
but at the cost of the low neutron detection efficiency 
and therefore much lower statistics. 



B. Measurement of p + 7r + coincidences 

We accept events with two charged reaction particles 
(p and 7r + ) in coincidence. They must show separate 
tracks in the wire chambers WC1, WC2, trigger sepa- 
rate sections of the E detectors, and at least one section 
of the F detector, but not the veto scintillator (V). The 
trajectories of the protons and pions are deduced from 
the wire chamber position readings. Their angular res- 
olution was limited primarily by multiple scattering in 
the 1.5 mm thick F detector and in the 0.18 mm thick 
stainless steel exit foil. Approximate angular resolutions 



(in the lab system) are a — 0.5° for protons and 1.0° for 
pions. This resolution was fully sufficient for the angular 
variations expected. 

The good intrinsic angular resolution of the wire cham- 
bers was used to check the consistency of the pion and 
neutron position readouts by tracing elastic protons to 
the hodoscope bars and comparing predicted and ob- 
served position readings. It was found that from run to 
run the beam axis and the detector symmetry axis could 
differ slightly in direction and also in their relative x and 
y coordinates at z=0 (the target center). We could also 
cross check the nominal z separation of the wire chambers 
since the separation and location of the hodoscope bars 
was fixed and well known. Small corrections of 1-3 mm 
had to be applied in software to the detector positions. 
After such corrections the remaining systematic angular 
error of the measured polar angles is about 0.04°. 

Charged particles that do not stop in or before the K 
detector trigger the V detector and are tagged as likely 
elastic events and generally vetoed. At 400 MeV we reach 
the design limit of the charged particle detectors and the 
veto detector begins to see (and reject) the most energetic 
pions at small angles. In deducing the energy spectra ac- 
count was taken of the differing non-linearity of light pro- 
duction for protons and pions by the plastic scintillators, 
as well as of energy losses in the exit foil, the F detector, 
air and other materials between the scintillators. After 
calibrations of all detector segments the detector stack 
provided an energy resolution for typical reaction pro- 
tons and pions of about AE/E = .09 (FWHM). 

The missing mass spectra contain a background contin- 
uum (see Fig. [l|) , which at higher beam energies stretches 
slightly beyond the missing mass peak. 
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FIG. 4. Neutron missing mass reconstructed from mea- 
sured pion and proton momenta. The background spectrum 
shown (dark area) results if the atomic hydrogen in the target 
cell is replaced by Nitrogen gas. 

The trajectory traceback indicates that this back- 
ground is primarily caused by beam halo hitting the Al 
and teflon components of the target cell. Without target 
gas (and the beam heating normally produced by it) al- 
most no background is seen. In order to obtain a realistic 
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background shape near and below the missing mass peak 
the target cell was filled with N2. This gas will produce 
some background of its own, but just as importantly it 
heats the circulating beam (as the Hydrogen gas would) 
and reproduces the ordinary beam halo. We found that 
the "iVa spectra" seen after the common software cuts 
looked identical to the background "tail" in the Hydro- 
gen missing mass spectra. Therefore, N2 spectra were 
measured with good statistics, and their shape was later 
used to correct for background under the missing mass 
peak. Our statistically most accurate measurements were 
obtained in the p+7r mode, i.e., by observing pions and 
protons in coincidence. 

C. Measurement of p+n coincidences 

Reaction neutrons in coincidence with protons were 
detected in a large hodoscope consisting of 16 long plas- 
tic scintillator bars. The bars were placed symmetrically 
about the beam direction in a plane defined by z=1.48 
m. They were 15 cm deep and mounted so that their 
dimension in the y and x directions were 120 cm and 5 
cm, respectively (see Fig. 1). The position in the y di- 
rection was determined from the differing arrival times 
of the scintillator light pulses read out by the top and 
bottom photomultiplicrs. The y-position resolution was 
a « 1.7 cm. At 325 MeV the geometric acceptance for p 
+ n detection is comparable to that for the pp — > ppir 
branch; however, the achievable event detection rate is 
much smaller because of the low neutron detection effi- 
ciency. The neutron pulse height threshold was set as low 
as practical and corresponds to 5 MeV electrons for all 
bars. At this threshold a 15 cm thick plastic scintillator 
averages a neutron detection efficiency of about 0.17 for 
the neutron energies of this experiment pTij] . 

A thicker neutron detector would be more efficient, but 
along with technical problems it would produce a cor- 
respondingly poorer time of flight resolution, since the 
length of the available flight path was limited to 1.5 m. 
In this experiment an additional reduction of the neutron 
detection efficiency arose because the E and K proton de- 
tectors are located in front of the neutron hodoscope and 
represent a 26 cm thick (polystyrene) absorber for the 
reaction neutrons. Resulting neutron losses in this "ab- 
sorber" range from 30% for the highest energy neutrons 
to about 90% for those at the very lowest energies. As 
a consequence the energy averaged effective neutron de- 
tection efficiency was reduced to a value of about 0.07. 
Since neutron energies are measured and neutron reac- 
tion cross sections are known pSf ], corrections for energy 
dependent efficiency losses can and have been made, but 
the loss in counting rate seriously limited the statistics 
obtained. 

The neutron energy was measured by neutron time of 
flight. In applying this method we use the correlated pro- 
ton trigger from pp — > pnn + in the F detector. Since the 



proton arrival at the F detector is delayed one has to use 
a two-step process: First, the trigger time difference (F 
detector time minus hodoscope mean time) is measured. 
Next the timing must be corrected for the proton flight 
time to the F detector, since the F detector is triggered 
by the proton after it has traveled about 30 cm before 
reaching the F detector. This correction is based on the 
measured proton energy and reconstructed track length. 
Neutron flight times (TOF) range from 5 to 12 ns. 

The dominant contribution to the TOF resolution 
comes from the 15 cm bar thickness, which constitutes 
10% of the flight path and cannot be overcome with 
the available detectors. Smaller contributions come from 
the intrinsic timing resolution of the hodoscope (0.4 ns 
FWHM) and the F detector (0.5 ns FWHM, after ampli- 
tude walk correction). We note that the raw time reso- 
lution of the F detector is worse than the figure quoted 
above because of trigger walk in the electronics and be- 
cause of the light loss and travel delay of light from parts 
of the large four-section F detector more distant from 
the photomultipliers. A substantial improvement was 
achieved by employing a pulse height compensation func- 
tion. Overall, we see a neutron time of flight resolution 
with AT/T fn 0.1. Therefore, the missing mass (MM) 
peak for tt + from p+n detection is not as sharp as for the 
corresponding neutron missing mass derived from p+7r + 
events. 



III. ANALYSIS 
A. Monte Carlo Simulations 

Our Monte Carlo (MC) simulations of the experiment 
used the event generator GENBOD of the CERN library. 
The simulation was used to determine various limiting 
effects of the apparatus and to derive corresponding cor- 
rections. The code contained the detailed geometry of 
the detector systems and the density distribution of the 
gas target. In the MC simulation we took into account 
the loss of energy of the charged particles before entering 
the detectors, detector resolutions, charged particle mul- 
tiple scattering, pion decay in flight, energy-dependent 
neutron detection efficiency and the probability of nu- 
clear reactions of the reaction neutrons in the E and 
K detectors. In the MC simulation we have used a pn 
final-state interaction (FSI) based on the Watson-Migdal 
theory, and the equations were derived following Morton 
p9| . We found that at the lower energies the FSI has a 
large effect on the overall coincidence acceptance. The 
simulation also provided a guide to the expected energy 
and angular distribution of the reaction products. 

Pion counting losses caused by the limited detector 
depth are not large enough to be detectable in the shape 
of spectra; however, the Monte Carlo simulation shows 
that they must be considered. At 400 MeV the loss for 
pions is 14% because this fraction of the forward pions is 
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too energetic to stop in the K detector. Only about 0.2% 
of the reaction protons penetrate past the K scintillators 
and are vetoed. The loss of high energy pions at small 
lab angles may create a small distortion of the 400 MeV 
p+7r data. Corrections to the 400 MeV spectra were not 
made since they would have to be very model dependent. 
We note parenthetically that the 400 MeV data from p+n 
coincidences do not have this systematic error, but within 
statistics they agree with overlapping p + 7r + results. No 
"veto" losses are seen at 375 MeV or below. 

The finite size of the individual detector segments pro- 
duces some counting losses since two sections have to 
trigger for acceptable events. However, systematic effects 
for the polarization observables are unlikely since the pro- 
tons have no strong <f> correlations with the pions. The 
segmentation used leads to a loss of about 7% in count- 
ing statistics for the p+7r + branch. There is no such loss 
for p+n detection. The charged particle detectors cover 
polar angles between 5° and 40° in the laboratory frame. 
Hence a large number of pions miss the detector. The 
total p + 7r + coincidence acceptance ranges from 21% at 
325 MeV to 15% at 400 MeV. 
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FIG. 5. Monte Carlo simulation for (a) the p+n and (b) the 
p+7r + acceptance, in the center of mass system. The partial 
acceptance for pions seen in the pn diagram at 0^ > 70° re- 
sults from the dominating forward boost for low energy pions. 
The cutoff at cos 9^ ~ —0.5 is caused by detector thresholds 
for the lowest ejectile energies. 

For p+n detection the MC simulation shows that the 
acceptance is symmetric about 90° although not quite 
isotropic. (See Fig. ^a.) Acceptance losses for p+n coin- 
cidences attributable to the detector geometry alone are 
of the order of 25%. The major cause is the central hole 
in the proton detectors. After all geometric acceptance 
losses and detector inefficiencies for neutron detection are 
taken into account the computed overall detection effi- 
ciency for pn coincidence events is 3.5%. It is seen in 
Fig. |^a that the angular variations of the coincidence ef- 
ficiency for the reconstructed pion are small. This is so 
despite the fact that we cannot detect protons at angles 
< 5° and neutrons at angles < 2.5° and have reduced 
coverage by the hodoscope of some azimuthal angles for 
large neutron polar angles. The Monte Carlo acceptance 
curves for p+n detection suggest that within the sta- 
tistical accuracy of the experiment the spin-correlation 



parameters integrated over 9^ and cV would need no sig- 
nificant correction. Fig. [|b shows the Monte Carlo sim- 
ulation for pn + acceptance as a function of cos 9 V . For 
pir + coincidences the apparatus acceptance is only useful 
for 9 V < 70°. Therefore, the integrated spin correlation 
coefficients will be deduced from the combined sets of the 
p + 7r + and p+n coincidences. 



B. Analysis of p+7r + coincidences 

The energies of the charged particles are measured by 
the plastic scintillator systems E and K. The calculated 
momenta of the unobserved particles strongly depend on 
the energies of the detected ejectiles, so considerable at- 
tention was given to a careful energy calibration of all 
detectors. The complex geometry of the segmented plas- 
tic detectors required corrections for light collection that 
primarily were derived from the observation of elastically 
scattered protons. An xy-position correction factor was 
applied to account for this dependence. A second pulse 
height correction factor was applied to compensate for 
a variation of phototube gains with the orientation of 
the magnetic guide field for the target polarization. For 



details see 22 



The corrected pulse heights L were converted into the 
deposited energy E using 



E = L + kiVL + ka 



(3.1) 



The nonlinear term corrects for light quenching in plastic 
scintillators, ko and k\ are calibration constants. L is the 
sum of the light pulse from the E and the K detectors 
in MeV and is given by L — ci(En ght + c 2 Ku g ht + c 3 ). 
The constants C\, C2, C3 are gain matching constants, and 
Eught and Ku g ht correspond to the observed light pulses 
in the E and the K detectors respectively. The constant 
C3 corrects for small energy losses in the material between 
the E and K detectors. It is small and set equal to zero 
when there is no K trigger. 

The total kinetic energy of the charged particle was cal- 
culated by also taking account of the energy lost by the 
charged particle on its way to the E-detector. The cali- 
bration constants were fine tuned by utilizing kinemati- 
cal relations. We required that the missing mass centroid 
was at its predicted value and that the angular distribu- 
tion of the pions from the simultaneous measurement of 
the reaction pp — > ppir was symmetric in the cm. sys- 
tem about 9^ — 90°. This symmetry was sensitive to the 
relative size of the calibration coefficients. However, the 
variation of the deduced spin correlation coefficients un- 
der different reasonable combinations of the calibration 
constants was small and less than the statistical errors. 
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FIG. 6. Detected pions at 325 MeV as a function of Q^. 
Only events with cos Q v > 0.4 were used for the analysis. 

Fig. ^ shows the directly observed 7t + differential cross 
sections plotted against cos 9^ in the cm. coordinate 
system. We note that there are almost no counts for 
pion back angles — 1 < cos 8^ < , as expected from the 
apparatus acceptance (compare Fig. ||b). 
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FIG. 7. Distributions of the calculated missing mass m x 
for p + 7v + detection at 325 MeV bombarding energy, for the 
four combinations of vertical beam and target polarization. A 
sharp peak (« 3.5 MeV/c 2 FWHM) is seen at 939.6 MeV/c 2 , 
the neutron rest mass. The shaded region indicates the back- 
ground distribution. 

Fig. ^ shows missing mass spectra seen at 325 MeV for 
four combinations of vertical beam and target polariza- 
tion at equal integrated luminosity. The polarization ob- 
servables are obtained from the ratio of "yields" for differ- 
ent spin orientations. The yields to be used are the inte- 
grated counts inside the missing mass gates minus back- 
ground. In order to estimate the error from uncertainties 



in the background we varied the background subtraction 
by ±25%. The effect on the final results was smaller than 
the statistical error. At 325 MeV the off-line resolution 
of the neutron missing mass peak was a = 1.4 MeV/c 2 . 
Even before software cuts and background correction it 
is apparent from Fig. that different spin combinations 
produce very different yields. 

It turns out that the decay in flight of pions plays a 
negligible role for these data. It will appreciably affect 
only the (undetected) backward scattered pions as these 
have much lower lab energies than the forward pions. For 
the final analysis we selected the pp — > pnir + events of 
interest by using a gate of 30 MeV or wider over the rele- 
vant missing mass peak. Gates as narrow as 10 MeV did 
not produce systematic changes, neither did they mea- 
surably reduce background induced errors. However, the 
narrower gates lead to some loss of statistics. 



C. Analysis of p+n coincidences 

This detection channel has the advantage that the ac- 
ceptance for the detection of p+n coincidences has little 
angular variation. So 0^ and ^-dependent acceptance 
corrections generally can be ignored. Therefore the p+n 
coincidences importantly complement the p+7r + channel. 
Reliance on the p+n angular distributions at large an- 
gles leads to larger statistical error bars relative to p+7r + 
(forward) region. However, the combination of the two 
detection modes provides data for the full angular range 
and so keeps the integrated spin correlation coefficients 
model independent. 

In the p+n analysis we first analyzed only those events 
where all three reaction particles (p, n and 7r + ) were de- 
tected (the triple coincidence). Next we evaluated the 
case where the pions missed the E detector, but a pro- 
ton and a neutron were detected (double coincidence). 
The energy of reaction protons was determined using 
the calibration constants described above. The energy 
of the neutrons was determined by measuring their time 
of flight (TOF) to the hodoscope. The MC simulation 
showed that although F was always triggered by pro- 
tons for a p+n double coincidence, in the case of a pmr + 
triple coincidence it was triggered by the faster pions. 
Therefore, depending on the event class, we corrected 
the neutron TOF by adding the time it takes either for 
the coincident proton or the pion to reach the F detector. 
A calculated offset was added to the timing signal of each 
hodoscope bar in order to make the timing information 
independent of the bar electronics. This correction was 
obtained by calibrating the timing circuits with elastic 
proton scattering. 

For 6 W} iab < 40° we observe pmr + triple coincidences, 
which are practically free of background. (The absence 
of accepted events from the N2 gas target showed that 
the triple pnn + hardware coincidence under standard 
software conditions eliminates all background from the 
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target wall and target impurities.) These events proved 
very valuable in assessing the correct shape of the miss- 
ing mass peaks in p+n and p+7r events. If a missing 
mass spectrum for triple coincidences is calculated based 
on the pion and proton momenta the (neutron) missing 
mass spectrum shows a very sharp peak as in Fig. 4, 
but there is no "background tail" at all. The triple co- 
incidence spectrum confirms the background subtraction 
shown in Figs. ^ and 

If the same triple coincidence events are used to calcu- 
late the (pion) missing mass by using the proton and neu- 
tron momenta (i.e., ignoring the simultaneously known 
pion momenta) we obtain the spectrum shown in Fig. |^. 
This spectrum can be used as a standard for the missing 
mass that p+n (double coincidence) events would have 
in the absence of background. 



2b0 1 




Missing mass (MeV/c 2 ) 

FIG. 8. Missing mass spectrum for pnn + triple coinci- 
dences, based on the measured neutron and proton energies. 
This spectrum contains no background from any competing 
reaction. The missing mass tail here is a consequence of some 
inaccurately measured neutron momenta. See text. 

The MM distribution peaks at the true pion mass of 
139.6 MeV, but there also is a "tail" over a wide range 
of the missing mass spectrum, which is not background 
related. We conclude that the counts in the MM tail of 
Fig. |8] represent genuine pnir + events from the Hydro- 
gen target, albeit events with poorly determined neutron 
momenta. We estimate that up to 20% of the p+n coin- 
cidences contain neutron observables that are distorted 
by interactions of neutrons with the K or E detectors. 
E.g., neutrons can undergo small angle elastic and in- 
elastic scattering, but still reach the hodoscope. This 
would lead to incorrect readings for polar and azimuthal 
neutron angles and hence to an incorrect missing mass 
calculation. 

Such events with poorly determined missing masses 
were excluded from further analysis. For all p+n events 
we reduce genuine background and avoid analyzing mea- 
surably distorted p+n events by using a missing mass 
gate from 100 to 160 MeV. 




FIG. 9. The ir + missing mass (m x ) spectrum at 375 MeV 
calculated from the measured neutron and proton momenta. 
The deduced background is shown by the lower distribution. 
For the analysis events with 100 < m x < 160 were accepted. 

Using the triple coincidence MM spectrum as a stan- 
dard, the background under the missing mass peak for 
two-particle p+n coincidences was deduced by adding a 
fraction of the measured unstructured V2 background 
continuum to the "standard" MM spectrum until the ob- 
served p+n MM spectrum shape was reproduced. The 
tail in the latter is flatter and more pronounced because 
of actual background contributions. To estimate the er- 
ror in this procedure we varied the match until it became 
unrealistic (±15%). A typical missing mass spectrum for 
p+n (double coincidence) detection is shown in Fig. ^. In 
the final result the uncertainty from this background sub- 
traction was about half as large as the statistical error. 
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FIG. 10. Energy and angular distributions for p+n coinci- 
dences at 375 MeV, compared with Monte Carlo projections 
(solid lines) in the laboratory system. Agreement is expected 
for the kinematic limits. However, the distributions may dif- 
fer because the 1^=0 assumption for the MC simulation is 
an oversimplification at and above 325 MeV. 

At 375 MeV the resolution of the pion MM peak was 
cr = 9 MeV/c 2 . Pion angular and energy distributions 
from p+n detection were computed using only events in- 
side this missing mass gate. Some resulting distributions 
are compared with Monte Carlo projections for the lab- 
oratory coordinate system in Fig. |l0|. 
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The end points of these distributions agree well with 
the kinematics of the experiment as they must. The solid 
curves represent pure 1^ = MC calculations. Although 
In = makes the major contribution, this MC assump- 
tion produces oversimplified energy and angular distri- 
butions. Nevertheless, the simulated distributions agree 
reasonably well with the data. 

Fig. [ll] shows the deduced pion angular distribution in 
the center of mass system. The reconstructed tt + distri- 
bution is plotted against cos 9 n in the center of mass (cor- 
rected for background and for the slightly non-uniform 
acceptance shown in Fig. ||a). As expected, it is non- 
isotropic and symmetric about 9^ — 90° within statistical 
errors. 
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FIG. 11. The relative tv + production cross section 
<To (cos 6 n ) at 375 MeV as deduced from p+n coincidences. 



IV. POLARIZATION OBSERVABLES 



A. Formalism for spin correlation coefficients 



The meaning of the symb ols A g used for polarization 

In terms of the "Carte- 



observables is defined by Eq. 4.1 



sian polarization observables" the spin-dependent cross 
section is written as 



<t(£,P,Q) = o-oK) 



i + J2 PMH) + £ QMS) + ]T FiQjA^) 



(4.1) 



where £ stands for the pion coordinates 9^,^^, the energy 
defining pion momentum^, and the proton coordinates 
9 P and <p p . The unpolarized cross section is o"o(£), and 
the polarization of the beam and the target is denoted 
by the vectors P = (P x ,P y ,P z ) and Q = (Q x ,Q y ,Q z ). 
The subscripts i and j stand for x, y or z and the sums 
extend over all possibilities. The resulting 15 polariza- 
tion observables include the beam analyzing powers Ajo, 
the target analyzing powers Aoj and the spin correlation 
coefficients Ay. 

The partial wave analysis for jyp pnir + is similar to 
that for pp — > ppir in terms of transition amplitudes. 







However, the pp — > pmr + transitions have isoscaler as 
well as isovector components. The different isospins in 
PP ~ * pnir + modify the selection rules for the reaction 
and lead to polarization observables that are different. 
The general relations between reaction amplitudes and 
angular distributions, however, remain almost identical. 
The applicable partial wave formalism has been discussed 
in detail in p6|| . We use the same notation as in ]26f| 
and reiterate some relevant definitions and theoretical 
relations below. Several names are in use for polarization 
observables. Their meaning is as defined below: 

A s (0 = A xx (0 + Ay V (0 (4.2a) 
A A (0 = A xx (£) - A yy (Z) (4.2b) 
Mt) = A xy (0~A yx (0 (4.2c) 

For identical particles in the entrance channel there are 
seven independent polarization observables: 

A y0 (0, A s (0, A zz (0, A xz (0, A A (0, A s (0 , A z0 (0. (4.3) 

This paper addresses the first five observables of this 
set. The remaining two, ^4s(£) and A z q(^), can be non- 
zero only for non-coplanar final states. In the following 
we will integrate over the angles of the nucleon, and thus 
these two observables vanish if parity is conserved. 

It is common to display the bombarding energy depen- 
dence of the observables in terms of the dimcnsionlcss 
parameter eta (77), which is defined as 

V = P-K,max/m n + , (4.4) 

The term "near threshold" is meant to include the energy 
region with 77 < 1, i.e., below 400 MeV. Setting c = % = 1, 
the maximum value of the 7r + momentum is found from: 

PTT,max = -^j=\J[s - (m n + m, p + m^+) 2 ][s - (m n + m p - m Tr +) 2 )] , (4.5) 



where y/s is the total center-of-mass energy, and m p , m n , 
and m^Jr are the masses of the proton, the neutron, and 
the pion, respectively. (We explicitly labeled the pion 
as 7r + to emphasize that the tt + and ir° mass difference 
matters here.) 

Below we quote some useful relations between inte- 
grated spin correlation coefficients and some directly ob- 
servable spin dependent cross sections. For two colliding 
spin 1/2 particles, one can define three total cross sec- 
tions, two of which depend on the spin. The total cross 
sections are related to the observables above by 

a tot = J aQiQdQpdQndpTr (4.6a) 
Aa T = - J a (£)Av(Z)dn p dn„dp v (4.6b) 
Aa L = -2 J a ({)A zz (Odn p dtt„dp n (4.6c) 

Here dH, = d cos 9 d<p, and the integration extends over 
< < 7r, and all pion momenta. Aa^/atd and 
Aax/crtot can have values between —2 and +2. 
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The integrated spin correlation coefficients are defined 
as: 



A s = [ 


1 (To {^)As(^)dQpdfl„dp^]/atot 


(4.7a) 


A zz = [ 


/ (To(g)A Z g(g)dQpdSlndpn]/cTtot 


(4.7b) 




1 cr (9^)Aa (9, ) sin9 7T d9,} / a tot 


(4.7c) 


A-xz ~ [ 


1 a (9,)A xz (9,)sin9,d9,]/(j tot 


(4.7d) 


AyQ = [ 


1 cr (9, ) A y o (9 n ) sinS, rf^] / a to t 


(4.7e) 



We note that A% and A& differ by a scale factor from 
Actt and Ao^ . These quantities in principle can be mea- 
sured directly, although in this study they are derived 
from integration over A^ (cos9,) and A zz (cos9,). The re- 
maining three integrals must be defined differently. Here 
the spin correlations A^ are taken at cj>n = 0. (They can- 
not be integrated over the variable (f>, since they would 
vanish, as will be seen below). 

Based on the dominance of Ss, Sp, Ps, and Pp tran- 
sitions, general symmetries and spin coupling rules p6| , 
the cross sections and spin correlation coefficients must 
have the general form: 

0"o (0 = aoo + &oo(3 cos 2 9^-1) + c (3 cos 2 9 P - 1) 

+do(3 cos 2 9, — 1)(3 cos 2 9 p — 1) + eo sin 29 p sin 29, cos Aip 
+/o sin 2 9 p sin 2 9 n cos 2Aip 

co(£)A/o(£) = [{a yQ + b yQ (3cos 2 9 p - 1) jsin^ + {c y0 + d y0 (3cos 2 9 p - 1)} sin 26^] cosip. 
+ [e y o + f y o cos9, -I- g y0 (3cos 2 9, - 1)] sin2# p cosip p 
+ [h y o sin^Tr + i y o s'm29,] sin 2 9 P cos(2ip p — ip„) 
+j y o sin2# p sin 2 9, cos^y^ — ip p ) 

M0^(0 = as + fe(3 cos 2 9,-1)+ c s (3 cos 2 9 p - 1) 

+ds(3cos 2 9 P — l)(3cos 2 9, — 1) + sm29 p sin 26*^ cos 
+/s sin 2 # p sin 2 9, cos 2 Ay? 

oo(£)4«(£) = a zz + 6^(3 cos 2 W - 1) + c Z2 (3cos 2 9 P - 1) 

+d zz (3 cos 2 P — 1)(3 cos 2 9, — 1) + e Z z sm2(9 p sin20,r cos Aip 
+f zz sin 2 6*p sin 2 9, cos 2A(^ 

o"o(0^a(0 = [a A + &a(3cos 2 6» p - 1)] sin 2 9, cos 2^ 
[ca + c?a(3 cos 2 6^ — 1)] sin 2 p cos 2<p p 
eA sin 2# p sin 29, cos((p p + ip,) 

Vo(QA xz (£) = [{a xz + b xz (3cos 2 9 p - l)}sin^ + {c xz + d xz (3cos 2 9 p - l)}sin2^ jr ] cos^ T 
+ [e xz + f xz cos9, + g xz (3cos 2 9, - 1)] sin20 p cos<p p 
+[h X z sm 9, + i xz sin 29,} sin 2 9 P cos(2</? p - <p w ) 
+j xz s'm29 p sin 2 9, cos(2tp, - <p p ) 
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Here we hctvG used ttic abbreviation A(p = epp — (p-n-- 

The 

Eqs. FS explicitly depend on the four angles 6 P , <p p , 9^ 
and ip„ . The energy dependent parameter p„ is contained 
in the coefficients. Statistics in this experiment are not 
sufficient to present double or higher differential cross 
sections. Therefore, we integrate over the angles of the 
proton and use energy and momentum conservation to 
eliminate all angles except V and This leads to a set 
of much simpler equations: 



*o(C) 


= a o + b 00 {3 cos 2 9^-1) 


(4.9a) 




= [a y o sin 9 V + c v q sin 29^\ cos ip„ 


(4.9b) 




= a s + 6 s (3cos 2 9 n - 1) 


(4.9c) 


o-o(CK*(C) 


= a zz + b zz (3cos 2 9 V - 1) 


(4.9d) 


*o(CMa(C) 


= a a sin 2 9-x cos 2<p-x 


(4.9e) 




= [a xz sin 9^ + c xz sin 29^] cos ip^ 


(4.9f) 



The symbol £ now represents the reduced set of vari- 
ables {ptt ,9^,4)^}. These equations display a simple and 
characteristic <f> n dependence of the different polarization 
observables and show the expected n dependence. The 
coefficie nts a v , b V) . . . for the set 4.9 correspond to those 
in Eqs. 4.8. They are obtained by one- or two-parameter 
fits to the observed angular distributions, separately for 
each observable. 



B. Extraction of polarization observables 

The data analysis, as described in the previous sec- 
tions, identifies the reaction particles, assesses the back- 
ground for each spectrum and calculates the kinematic 
variables and spin dependent cross sections of the reac- 
tion products. It produces event files which contain kinc- 
matically complete information for all detected reaction 
particles. For each beam energy there are 12 such event 
files, one for each combination of beam and target spin. 
These yields are first corrected for the beam luminos- 
ity, which can vary for beam "spin up" and "spin down" 
subcycles, and for the background measurement. The 
background correction was made for each selected 9^ an- 
gle bin individually. The ratios Ri of yields for different 
spin combinations, integrated over a chosen 9„ range, are 
then analyzed as a function of <p w , because the allowed 
6jr dependence can be predicted from spin coupling rules 
||2lf . For this energy range, only final states with pn or 
pion angular momentum of and 1 are expected to be 
significant. In a previous measurement of pp — ► dn + at 
400 MeV it was found that any l n = 2 contribution is 
very small f30|| . This allows us to consider only transi- 
tions to Ss, Sp, Ps and Pp final states in the analysis. Sd 
and Ds transitions would affect the energy dependence 
of the coefficients only and so are very difficult to sep- 
arate from Pp transitions |2q|. They will be ignored in 



this analysis. We then have explicit pr edic tions for the 
expected 9 and <fi dependences from eq. 4.9 . 

The combination of p+7r + and p+n measurements pro- 
vides model-independent values for the polarization ob- 
servables for all polar and azimuthal angles of the pion. 
The low neutron detection efficiency and the resulting 
low statistical accuracy of the p+n data make it advis- 
able to display the combined data using some theoretical 
guidance. As shown below, Ay,{9^), A^{9„), and A zz (9 7r ) 
must be symmetric about 9 V = 90° for the transitions 
considered. So a good analysis in terms of the pion co- 
ordinates does not require the (redundant) data at large 
polar angles. This simplification and the fact that all 
published theoretical predictions have been presented in 
terms of the pion coordinates makes these coordinates 
our preferred system for the analysis. 

The microscopic relations between the coefficients and 
the transition amplitudes can be derived from the partial- 
wave expansion described in the appendix of [|26| , but 
they are complicated. Moreover, the number of indi- 
vidual pp — > pnir + amplitudes contributing above 350 
MeV has become too large (19 rather than the 12 for 
pp — > ppir ) , since isospin 1 and are present in the final 
state. They could not be deduced individually from the 
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FIG. 12. The yield related ratios (Ri - l)/{Ri + 1) as a 
function of the pion azimuthal angle 4>n for data integrated 
over all other coordinates. The specific beam and target spin 
combinations selected are listed on the left. For multiple ar- 
rows the orientation of the beam spin is shown first. Longi- 
tudinal polarization is indicated by the symbols o (opposite) 
and a (along beam direction). The solid curves represent a 
least-square fit using the expected theoretical (j>n dependence. 

When ca lcul ating the value of a polarization observable 
from Eqs. 4.8 or 4.£, one evaluates the ratio Ay(£) = 
(Tq (£) Aij (£) /cr (£) , so the overall normalization of all 
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terms in these equations cancels. As seen from Eqs. 4.9 
the yield ratios Ri(4>-n) could either be constant or have 
a <p v or 24>tt dependence. This is borne out by the data 
(compare Fig. |l2| ). 

The polarization observables were deduced by evalu- 
ating the observed 4> v dependences of the ratios Ri for 
selected beam and target spin combinations. This evalu- 
ation is complex when longitudinal as well as transverse 
beam polarizations are present at the same time. There- 
fore, the devolution process uses the computerized fitting 
routine BMW pl| , which was written for this purpose. 

Fig. |l2| shows the (f>^ dependence of six spin dependent 
yield ratios. The data for the beam (first arrow) and 
target spin combinations indicated have been integrated 
over all coordinates other than the coordinate <f> w . The 
cur ves are fits using one to three components of the Eqs. 
4.9 . The first three rows present different ways to extract 
the analyzing power A y (4>tt)- The lower three rows con- 



tain information on A-^ = A x 



A,.,,, A a — A„ 



and contributions from A zz and A xz . For some ratios 
the statistics are marginal, and one cannot exclude the 
potential presence of I components higher than included 
in the analysis, but the cj)^ dependent fits show that the 
inclusion of Ss, Sp, Ps, and Pp transitions is sufficient to 
reproduce the data within experimental errors. 

For the simultaneous detection of neutrons and pro- 
tons, our data sample the full range for 9^, although 
with low statistics. We combine the p+7r and p+n data 
sets to obtain optimal spin correlation coefficients for the 
full angular region. We avoid difficulties generated by 
the non-uniform detector acceptances in 9^ by evaluat- 
ing the p+n and p + n relations A y (cos 9^), which are 
ratios of cross sections at a given angle. (Our detection 
efficiency does not depend on spin, and the detector ac- 
ceptances cancel out for Aij(cos 9^).) The combined p+n 
and p+7r + sets yield complete angular distributions with 
their best statistics at forward angles. The unpolarized 
angular distribution ao(9 v ) was obtained to sufficient ac- 
curacy from the p+n branch. The angular distributions 
can now be integrated. To best account for experimental 
errors we have chosen to integrate Eqs. 4.9 directly after 
the fitting coefficients are deduced. 

Some of the polarization ratios measured are not inde- 
pendent, as the first three rows in Fig. 12 show. The re- 
action has additional redundancies. If parity is conserved 
and if we have identical particles in the entrance channel 
this redundancy can give us back-angle information for 
Ay even though our detectors only cover forward polar 
angles for the direct detection of pions. The correlation 
we use repeatedly is 



Aij(9 p , 4> p , 6^, (f>„) = Aji(ir - 



+ 7T, 7T 



(4.10) 



This relation holds for i ^ j and also for i=j. E.g., since 
both A xz and A zx are measured at forward angles, we 
will obtain the back angle information for A xz from the 
A zx measurement at forward angles. The polarization 



observables A xz and A zx are not symmetric about 9 n 
90°, so this redundancy becomes very useful. 



V. RESULTS 



A. Polarization observables 



It follows from Eq. 4.1C that the observables A-^, A&, 
and A zz are symmetric about 9^ — 90° (cos 9 n = 0) . 
Within statistical errors the experimental data agree with 
this expectation. In Fig. 13 we have reduced the scatter 
from the low statistics of the p+n coincidences by com- 
bining the corresponding data for forward and back polar 
angles. The data for cos^ > 0.5 are dominantly deter- 
mined by events from p + n coincidences. In agreement 
with theoretical expectations there is only a slow depen- 
dence on the polar angle, so the lack of good statistics 
near 9 = 90° does not impede comparison with theory or 
the extraction of good values for the integrated polariza- 
tion observables. 

Fig. [l4| shows results for A y (9 7r ) and A xz (9 7r ) for the 
full angular range so potential asymmetries can be seen. 
The statistically most accurate data were obtained for 
375 MeV. Here and at 400 MeV the Jiilich model is at 
odds with the data. The fit with Eq. |4.9| (solid lines) 
does much better. Still, a close inspection of the fits 
shows some small, but statistically significant differences 
between the partial wave curve and the data at very small 
and very large angles. We also see from Table I that the 
X 2 value for the A y fit has become large. The A y data 
suggest that higher partial waves enter at 375 MeV, but 
the experimental uncertainties discourage the extraction 
of relatively small contributio ns. 

The fits obtained with Eq. 4.9 are good (i.e., \ 2 ~ 1 
for all cur ves except for A y at 375 and 400 MeV). There- 
fore Eqs. 4.9 together with the coefficients of Table II 
can be used to represent the new data. The coefficients 
in these equations are bilinear sums of the reaction am- 
plitudes. Their experimental values are given in Table 
II. This set is also used to obtain the integrated spin 
correlation coefficients. Integration of the angular distri- 
butions shown above produces the spin correlation coef- 
ficients in Cartesian coordinates. These coefficients were 
the original objective of this experiment. They are now 
known with good statistical accuracy and are given in 
Table III. A comparison of these integrated polarization 
observables as a function of beam energy with predictions 
of the Jiilich model is shown in Fig. [15| 

For completeness we note here that our attempt to 
extract non-coplanar angular distributions for A z q pro- 
duced results (not shown) that were consistent with zero. 
This is different from the pp — > ppir results, but has no 
obvious explanation. 
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FIG. 13. 8tt dependence for the polarization observables As = A xx + A yy , Aa = A xx — A yy , and A zz , in cartesian units. 
Data for the range 0.5 < cos 8 < 1 are primarily determined by the p+7r coincidences. The remaining points come from the 
p+n coincidences. The error bars include all random errors as well as estimated uncertainties from background subtraction. 
The dashed lines are Jiilich model predictions by Hanhart et al.[22]. The solid lines show fits with the Eqs. ft.Sj. 
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It is clear from Figures |l3|, [lj, and [T| that the dis- 
tributions based on the Julich model are in good agree- 
ment with the data at 325 MeV. However, above 77 = 0.7 
they produce ever larger \ 2 values when compared to the 
data. These disagreements become striking for and 
A y . The failures are most visible for A y , an observable 
sensitive to admixtures of higher partial waves. (More 
serious disagreements with this model have been seen for 
the isovector production in pp — » ppit p6[ . However, 
as discussed below, in this energy region isovector terms 
contribute less than 10% to the pp — ► pmr + cross section. 
The observed differences in pp — > pnir + grow well beyond 
this level.) 

Our partial wave analysis, which includes Ss, Sp, Ps, 
and Pp transitions, generally provides fits to the mea- 
sured angular distributions with % 2 (per degree of free- 
dom) values near 1. The exceptions are A y at 375 and 
at 400 MeV, where the cross sections are largest and the 
statistics are good. Some \ 2 values as large as 3.9 are 
found if only statistical errors are considered. 

The values for the product P*Q are known to good 
precision (see Table I), but errors for the beam (P) or 
target (Q) polarization individually are not negligible at 
the lower energies. Changes in P and Q affect only the 
analyzing powers A y (6). They could reduce or increase 
the asymmetry of the angular distributions. Typically, 
the uncertainties in P arc smaller than the statistical er- 
rors. 



B. Discussion and comparison with other work 

The statistical and fitting errors listed in Tables II 
and III include all known and estimated random errors. 
As explained above, all angles were measured simulta- 
neously, and systematic normalization errors for Aij are 
unlikely Based on the detector design and redundant 
measurements we expect that all systematic errors have 
remained small. In the center region (cos 6 w 0) the 
angular distributions show large statistical errors. How- 
ever, these data points do not materially affect the par- 
tial wave fits or the integrals. We note that our initial 
results reported in Q were subject to some model de- 
pendence that is absent here. Nevertheless, they are con- 
sistent with the final results presented here. Noticeable 
asymmetries around 90° have been seen for A y above 
350 MeV. In the framework of our partial wave analysis 
this asymmetry must be produced by Pp transitions. (At 
higher energies such asymmetries can also be produced 
by Ds and Sd transitions.) With the possible exception 
of the analyzing powers A y {&) at 375 and 400 MeV, the 
pion production data are well represented by the partial 
wave predictions based on the the assumption of Ss, Ps, 
Sp, and Pp transitions. The Julich model predictions 
and the data agree for A&. However, we see serious dis- 
agreements for As and A y as the beam energy increases. 
Ref. jll| includes more amplitudes than our analysis, but 



the calculations predict little asymmetry for A y (9). The 
differences for A y and the increasing divergence with en- 
ergy are also seen in Fig. 15. At this time there are no 
predictions available for A xz and A z0 . 
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FIG. 15. Energy dependence of the integrated spin corre- 
lation coefficients An, Aa, A zz A xz and A y o, and the peak 
analyzing power A ytm ax for pp — > pnn + . The diamond shape 
symbols represent measurements at lower energies and are 
taken from Ref. [18]. The solid lines are predictions of the 
Julich meson exchange model. (There is no prediction for 

A X z- ) 



C. Deduction of important partial waves 

The number of contributing partial waves grows 
rapidly with energy. If we restrict ourselves to Ss, Sp, 
Ps, and Pp contributions as above, the 19 individual am- 
plitudes listed in Table IV are needed for a detailed inter- 
pretation of the data. There are 12 isoscalar and 7 isovec- 
tor amplitudes. The experimental information available 
includes the three cross sections a tot, Aut, and A ox for 
pp — > pmr + that are related to these amplitudes. In addi- 



tion a recent pp — > ppir study provides the three relevant 
isovector cross sections cx' tot , Acr' T , and Ao' L (Ref. p3| , 
Table V). As long as isospin is a good quantum number 
these cross sections also give the isovector part of the 
pp — ► pnir + reaction if taken at the same 77. 
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So one has 6 new measurements for 19 variables. This 
necessitates some restriction of the further analysis. In 
a previous pp — > pnir + study jll| , closer to threshold, 
the partial wave space had been restricted to the lowest 
isoscalar amplitudes ao,<2i,<X2, an d to the lowest known 
isovector amplitudes. With this simplification and with 
reliance on the measured analyzing powers three ampli- 
tudes were deduced for r\ < 0.5. Some of these earlier 
results will be shown below. It will become apparent 
in comparison with our new data that the angular mo- 



mentum space considered in Ref. |18| is too small for 
7] > 0.3. For 0.3 < rj < 0.9 it becomes necessary to 
consider all Ss, Sp, Ps, and Pp contributions. In order 
to reduce the number of variables we use similarities in 
the spin algebra coefficients for the 19 amplitudes of in- 
terest. A suitable combination of the 19 partial cross 
sections into six groups allows us to find the Sp and Ps 
strengths separately to deduce the lowest Pp isoscalar 
partial cross section for the amplitude 0,3 directly, and 
to put a close upper limit on the Ss contributions. We 
will identify the isoscalar partial wave cross sections by 
cr(ao), cr(ai), <j{a,2), ■ ■ ■ and the isovector partial wave 
cross sections by <r(b ), a(b{), (7(62), ... as in Table 
V. Generally, a(ai) = Ci\a,i\ 2 , where the C, factor is a 
combination of factors like it and Clebsch-Gordan coef- 
ficients, which can differ from amplitude to amplitude. 
(Therefore the partial cross sections listed in Table V do 
not provide the magnitude of the corresponding ampli- 
tudes without further work.) The notation a{a\,a^o) 
implies that we could not separate the cross sections for 
01, the Ss component, from the Pp components 04 to 
clq. Hence cr(ai, <i4->6) = + c(a4) + (7(05) + a(ae). 

The partial cross section groups that could be isolated 



are given in Eq. 5.1 



Sp isoscalar terms : c(ao, 0-2) = r (Actl + 2Aot + 2cr tot — Aa' L — 2Ao' T — 2a' tot ) (5.1a) 

8 

Ps isovector terms : c(&i, 62) = - {Aa' L + HS.a' T + 2a' tot ) (5.1b) 

8 

Ss + Pp isoscalar terms : <r(ai, 04^) = t (— Ao^ + 2otot + Aa' L — 2a' tot ) (5.1c) 

Ss + Pp isovector terms : cr(b ,b 3 ) = ^ (Aa' L - 2Aa' T + 2a' tot ) (5. Id) 

8 

Pp isoscalar terms : 17(03) = r (Ao\l — 2Actt + 2a to t — Aa' L + 2Ao' T — 2a' tot ) (5.1e) 

8 

Pp isovector terms : c(&4^ii) = ^ (— Aa' L + 2cr' tot ) (5.1f ) 



Of these six equations, which hold for pp — > pnir + , 
three also hold for pp — * ppir . We note that Eq. 5.1b has 
been presented before. It is identical to Eq. 13 in p6| . 
The six equations now permit a calculation of partial 
wave cross sections to the specified groups of final states 



1G 



from the measured spin-dependent cross sections. The 
sum of these partial cross sections equals the total 7r + 
production cross section. Since the partial cross sections 
add incoherently the effect of higher lying weak ampli- 
tudes is minimized. This is an advantage over relying 
on analyzing powers, which are sensitive to even small 
admixtures. The amplitudes incl ude d in each group are 
indicated on the left side of Eqs. 5.1. 



In many experiments, including the present one, it 
is much easier to measure accurate cross section ra- 
tios than absolute cross sections. So our experimental 
pp — > pmr + quantities are given as a fraction of the total 
7r + production cross section a tot- The Eqs. |5.1| are easily 
re-written in terms of partial wave strengths by dividing 
both sides by a tot- For the p + tt + branch the values 
Aax/o'tot and Ac7£/<7 tot were calculated from eqs. 4.6 
and As(0), Aa(9) and <tq(6). In figures and tables we 
will generally use the ratios of partial wave cross sections 
to total cross sections. We refer to them as partial wave 
strengths. 
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FIG. 16. Partial wave strengths for six groups of ampli- 
tudes as function of r\. The isoscalar cross sections are con- 
nected by solid lines, the isovector ones by dashed lines. The 
contributing amplitudes, including the (small) Pp contribu- 
tions not resolved from the dominant Ss cross sections are 
indicated in the legend. The dash-dotted line represents the 
full Pp isovector strength contribution in pp — > pmv + , as de- 
rived from the results of Ref . [26] . 

To use them in this study the total pp — ► ppir and 
pp — ► pmr + cross sections were taken from the literature 
and interpolated for the present 77 values. We obtained 
the pp — ► ppir information needed from Ref. p6| and the 
pp — > pmr + total cross sections from and from Fig. 2 
in Ref. || . The accuracy of the total cross section ratios 
so obtained is not very high, but it will suffice here be- 
cause the isoscalar terms of interest are an order of mag- 



nitude larger than the isovector ter ms. T he partial cross 
section strengths derived with Eqs. |5.l| are displayed in 
Fig. 16 and listed in in Table V. 

The primed cross sections are the (pure isovector) cross 
sections measured for pp — > ppir , which are also more 
accurately given as fractional strengths. To work in 
terms of pp — > pmr + partial wave strengths the pp — > 
pp-K° strengths of Ref. |2(| have to be multiplied by the 
ratio of the pp — > ppir and pp — > pnir + unpolarized cross 
sections, taken at the same relevant 77 values. 

Fig. |l6| shows the change of partial wave strength with 
energy for Sp, Ps, and other groups. It is immediately 
apparent that for the energy region studied the leading 
isoscalar partial cross sections are an order of magni- 
tude larger than the isovector ones. It helps our discus- 
sion that lowest-lying Pp isoscalar partial wave strength 
Pp(aa) could be resolved. It is much smaller than the Ss 
and Sp strengths. So is the sum of all isovector cross sec- 
tions for b4 to bn. The Pp strengths attributable to 63 
can be assessed by comparing Pp(64_>n) from this work 
with the heavy dash-dotted curve derived from pfj for 
the full Pp isovector strength Pp(6 3 _ > n). 
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FIG. 17. Sums of isoscalar and isovector partial wave 
strengths as function of 77. The Sp+Ps sum is measured di- 
rectly. The points labeled "Ss max" represent a close upper 
limit to the sum of the Ss partial cross sections. Any cor- 
rection for the unresolved Pp amplitudes (a4,a5,aa, and 63) 
would lower the Ss curve (as indicated by the estimated er- 
rors). The admixtures can be expected to be smaller than 
Pp(as). The lower points show the documented Pp strengths 
only. The data at r) = 0.22 and 0.42 are from Ref. 18. 

Therefore, it seems reasonable to assume that the Pp 
contributions from 0,4, (15, a§ and 63, which could not be 
disentangled from the the Ss amplitudes, are also much 
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smaller than the Ss terms. On this basis we estimate 
that they make up no more than 5 to 10% of the "Pp 
entangled" Ss max cross section curve. For r\ = 0.9 the 
Sp(ao,et2) cross section has become dominant. As seen 
in Fig. 16, it is very much larger than the Ps isovec- 
tor contribution. It would be of interest to resolve the 
isoscalar component do, because it can be used to con- 
strain the strength of 3-body forces jljj. However, in 
this analysis do and the much larger amplitude 02 always 
appear together. The Ss fraction, including the unre- 
solved Pp contributions, has fallen to less than 0.3. This 
is consistent with the work at 420 MeV |3S|. 

In Fig. [l7| the data points give the summed Sp+Ps 
strengths, the upper limit for the summed Ss strengths, 
and a lower limit for the Pp strength. The heavy dashed 
curve shows the likely energy dependence of the actual 
Ss strength. The divergence of the old and new pp — > 
pmr + interpretation near rj « 0.45 serves as a reminder 
that a partial wave analysis is only model independent 
if it fully encompasses all contributing amplitudes. This 
apparently was no longer true for the 320 MeV data (77 = 
0.42) of Ref. 18. 

In this respect our present difficulty to perfectly repro- 
duce A y at 375 and 400 MeV in the Ss, Sp, Ps and Pp 
frame work (see Table II) should be taken as a warning. 
At these energies some higher partial waves may con- 
tribute enough so that they must be considered, at least 
for the analyzing powers. 

VI. CONCLUSIONS 

We have measured the spin correlation coefficients 
A'E=A X x~\~Ayy J A^—A X x~Ayy^ A zz , A xz , and Ay, as well 
as angular distributions for cr(0 w ) and the polarization 
observables Aij{0 n ) at energies from 325 to 400 MeV. At 
the lowest energies the results are in agreement with pre- 
diction of the Julich meson exchange model. The agree- 
ment deteriorates considerably at energies where Ss tran- 
sitions no longer dominate. At 375 and 400 MeV some 
physics aspects in pp — > pnir + apparently are missed 
by the model. This suspicion is supported by the even 
poorer agreement of the model with the pp — > ppir data 

H- 

The pp — > ppir and pp — > pnir + reactions are found 
to differ greatly in the relative importance of Sp, Ps, and 
Pp transitions. Sp strongly feeds the delta resonance in 
pp — > pnn + , but this transition is forbidden for pp — > 
ppir . By contrast Ps contributions in pp — ► pnir + are no 
larger than Pp contributions, as seen in Fig. 16. In pp — > 
pmr + the Ss and Sp isoscalar terms are most important 
while the Pp transitions just begin to contribute. For 
pp — > ppir Pp becomes dominant at rj = 0.7. 

The partial wave analysis was able to reproduce almost 
all polarization observables within experimental errors. 
This supports the postulated adequacy of considering 
only Ss, Sp, Ps, and Pp transitions in the near threshold 



region. However, this angular momentum space may not 
be adequate to explain details of analyzing powers, be- 
cause they can be affected by small admixtures of higher- 
lying transitions. Even in this limited space the number 
of individual partial waves for pp — ► pnir + at 400 MeV 
is too large to deduce all individual amplitudes. Some 
interesting s um rules for groupings of amplitudes were 
found (Eqs. |5.l| ), and the corresponding partial cross 
sections could be extracted. They show, e.g., that for 
rj < 1 Pp and Ps terms play a considerably smaller role 
in pp — > pnir + than in pp — ► ppir . 

Further progress may come from improved theoretical 
models that can accurately predict the new data at hand. 
It is interesting to note again that the Julich model does 
well at 325 MeV where Ss dominates, but it increasingly 
fails for pp — > pnir + (as well as for pp — > ppir ) as higher 
angular momenta become important. 
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TABLE I. Beam energies, integrated luminosities for the p + tv + measurements, and the products of beam and target 
polarization for runs a and b. (No p+n data were taken in run a. The p+n measurements began in the middle of run b and 
have correspondingly lower integrated luminosities.) 





Run a 


Run b 


Energy 


j Ldt 


PyQ 


j Ldt 


P X Q 


PyQ 


P Z Q 


(MeV) 


(nb- 1 ) 




(nb- 1 ) 








325.6 


2.163 


0.456 ±0.003 


3.0 


0.059 ±0.002 


0.333 ±0.002 


0.296 ±0.003 


350.5 


0.901 


0.342 ±0.004 


1.3 


0.053 ±0.003 


0.316 ±0.003 


0.267 ±0.005 


375.0 


3.024 


0.514 ±0.004 


4.1 


0.041 ±0.002 


0.333 ±0.002 


0.266 ±0.004 


400.0 


0.831 


0.526 ±0.006 


1.1 


0.039 ±0.004 


0.289 ±0.004 


0.203 ±0.008 



TABLE II. Coefficients for the fits with Eqs. 4.9 that reproduce the measured angular distributions of the polarization 
observables. The associated Legendre polynomials used for the fits are determined by selection rules for Ss, Sp, Ps, and Pp 
transitions. The unpolarized angular distribution 0o(aoo,boo) is given in arbitrary units by setting aoo = 1. The errors listed 
refer to the individual fitting coefficients. The \ 2 numbers give the overall quality of the fit to the data per degree of freedom. 
The fits are shown in Figs. 11, 13 and 14. 



Name 


par am. 
value 


325 MeV 
par am. 

error 


x 2 

data 


par am. 
value 


350 MeV 
param. 
error 


2 

X 

data 


param. 
value 


param. 
error 


375 MeV 

x 2 

data 


param. 
value 


param. 
error 


400 MeV 

x 2 

data 


aoo 
boo 


1 

0.168 


0.035 




1 

0.190 


0.040 




1 

0.199 


0.030 




1 

0.196 


0.045 




as 
b s 


-0.560 
-0.303 


0.052 
0.063 


0.5 


-0.810 
-0.478 


0.055 
0.067 


0.6 


-0.994 
-0.510 


0.015 
0.018 


2.8 


-1.070 
-0.439 


0.024 
0.029 


1.4 


(Ja 
b A 


-0.037 


0.091 


1.4 


-0.001 


0.097 


0.2 


0.084 


0.028 


1.3 


0.075 


0.045 


1.4 


a zz 
b zz 


0.120 
-0.188 


0.042 
0.054 


0.1 


-0.177 
-0.257 


0.047 
0.057 


0.8 


-0.310 
-0.233 


0.018 
0.021 


1.0 


-0.431 
-0.199 


0.037 
0.043 


3.1 


a y o 

CyQ 


-0.247 
0.007 


0.015 
0.013 


0.5 


-0.255 
0.050 


0.013 
0.010 


0.9 


-0.276 
0.044 


0.005 
0.005 


2.4 


-0.285 
-0.032 


0.008 
0.007 


3.9 


Oixz 
C-xz 


-0.051 
0.106 


0.042 
0.038 


0.7 


0.021 
0.047 


0.042 
0.036 


1.2 


0.053 
-0.040 


0.021 
0.020 


1.1 


-0.041 
-0.104 


0.041 
0.036 


1.2 



TABLE III. 


Beam energy, the r\ parameter, 


and the deduced intej 


grated spin correlation coefficients. 


The table gives the 


weighted average of all 


runs as shown in Fig. 15 










T(MeV) 


V 




Aa~ 


A zz 


Ay 


A xz 


325.6 


0.464 


-0.533 ±0.046 


-0.027 ±0.064 


0.148 ±0.041 


-0.209 ±0.011 


-0.043 ±0.044 


350.5 


0.623 


-0.761 ±0.046 


0.001 ±0.070 


-0.143 ±0.040 


-0.218 ±0.011 


0.018 ± 0.036 


375.0 


0.753 


-0.945 ±0.068 


0.062 ±0.036 


-0.283 ±0.016 


-0.237 ±0.004 


0.045 ±0.019 


400.0 


0.871 


-1.026 ±0.023 


0.056 ±0.034 


-0.414 ±0.032 


-0.244 ±0.011 


-0.035 ±0.036 
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TABLE IV. Angular momentum quantum numbers for the 
partial waves of the reaction pp — > pnir + . 



Type 




Label 


Os.Xli 9c 1 -4-1 1 .1 


Ss 


isoscalar 


ffli 


tt; 'j _ 

3 Pi ^ 3 Si,s 


Ss 


isovector 


bo 


3 P ^S ,s 


Sp 


isoscalar 


ao 


x 5 ^S^p 






0,2 


1 P>2 -> 3 Sl,p 


Ps 


isovector 


bi 


l S ^ 6 P ,s 






b 2 


1 D 2 ~^ 3 P 2 ,s 


Pp 


isoscalar 


a.3 


3 P ^ L Pi,p 






04 


3 Pi ~> 1 P 1 ,P 






ffls 


3 P 2 ^ 1 P 1 , P 






CS6 


3 F 2 ^Pi,p 


Pp 


isovector 


h 


3 P ^ d Pi,p 






bi 


3 P 2 -> 3 Pi,P 






bs 


3 P 2 ^ 3 P 2 ,P 






be 


3 P 2 ^ 3 Pi,P 






br 


3 P 2 -> 3 P 2 ,P 






b s 


3 Pi ^ 3 P ,p 








3 Pl -^ 3 P!,p 








3 Pi -^ 3 P 2 ,P 








3 F 3 -> 3 P 2 ,p 



TABLE V. Listing of the pp — > pnir + partial wave strengths for the groups of isoscalar and 
isovector amplitudes indicated. The 300 MeV results listed were taken from Ref. [18]. The 300 
MeV strengths not listed are assumed to be negligible. 



Isoscalars 
















E(MeV) 


V 


Sp(a0,a2) 


error 


Ss+Pp (al,a4-6) 


error 


Pp(a3) 


error 


300 


0.220 


0.088 


0.030 


0.740 


0.050 






325.6 


0.464 


0.342 


0.018 


0.570 


0.024 


0.006 


0.018 


350.5 


0.623 


0.469 


0.018 


0.421 


0.023 


0.052 


0.018 


375 


0.753 


0.541 


0.020 


0.342 


0.011 


0.050 


0.020 


400 


0.871 


0.579 


0.013 


0.262 


0.019 


0.063 


0.013 


Isosvectors 
















E(MeV) 


n 


Ps(bl,b2) 


error 


Ss+Pp (b0,b3) 


error 


Pp(b4-ll) 


error 


300 


0.220 






0.173 


0.022 






325.6 


0.464 


0.007 


0.002 


0.073 


0.007 


0.004 


0.003 


350.5 


0.623 


0.014 


0.003 


0.044 


0.005 


0.007 


0.002 


375 


0.753 


0.019 


0.002 


0.034 


0.003 


0.016 


0.002 


400 


0.871 


0.024 


0.003 


0.034 


0.004 


0.031 


0.004 
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